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ABSTRACT 


This  report  documents  two  computer  programs  which  determine  the 
wave  resistance  and  propulsive  performance  of  Small-Waterplane-Area 
Twin-Hull  (SWATH)  ships.  These  programs  are  SYNTHESIS  and  CLOSEFIT, 
and  they  use  as  input  data  the  moments  and  offsets  of  ship  geometry, 
respectively.  While  they  both  employ  linearized  ship  wave  theory  and 
thin-shlp  approximation  and  share  many  subroutines,  CLOSEFIT  is  a  more 
refined  program  capable  of  furnishing  superior  results.  On  the  other 
hand,  SYNTHESIS  is  simpler  and  includes  several  additional  subroutines 
which  can  determine  the  characteristics  and  performance  of  a  SWATH 
ship's  propulsive  system  at  both  design  and  off-design  conditions. 


ADMINISTRATIVE  INFORMATION 

This  investigation  was  authorized  under  a  direct  funded  block  from  the  Naval 
Material  Command  (NAVMAT  08T23)  under  Program  Element  62543N,  Task  Area  ZF43-421- 
001,  and  funded  through  the  High  Performance  Vehicles  Office  of  the  David  W.  Taylor 
Naval  Ship  Research  and  Development  Center  (DTNSRDC),  Work  Unit  1-1500-103. 
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BACKGROUND 


This  report  (Maintenance  Manual)  contains  documentation  for  two 
computer  programs,  CLOSEFIT  and  SYNTHESIS,  which  have  been  developed  at 
the  David  W.  Taylor  Naval  Ship  R&D  Center  (DTNSRDC)  to  predict  the 
effective  and  delivered  power  for  Small-Waterplane-Area  Twin-Hull  (SWATH) 
ships.  This  maintenance  manual  contains  documentation  for  each  subroutine 
in  the  programs,  definitions  of  all  of  the  variables  in  the  COMMON  BLOCKS, 
and  a  listing  of  the  computer  programs.  This  documentation  is  intended 
for  the  individual  who  must  modify  the  computer  programs,  and  must 
therefore  be  able  to  understand  the  functioning  of  the  programs,  and 
should  be  contrasted  with  the  Users  Manual,  which  serves  to  define  the 
input  to  and  output  from  the  programs.  The  ability  to  prepare  proper 
input,  and  to  use  the  output  does  not  require  an  understanding  of  the 
workings  of  the  program. 

SWATH  ships  are  catamarans  which  combine  a  thin  strut  with  a  submerged 
elongated  slender  body.  Each  demi-hull  of  the  SWATH  has  either  one 
(single)  or  two  (tandem)  struts  which  support  the  lower  hull.  Figures  1 
and  2  show  schematic  diagrams  illustrating  the  geometry  of  single  and 
t_  «aem  strut  SWATH  configurations,  respectively. 

Programs  CLOSEFIT  and  SYNTHESIS  have  a  similar  theoretical  foundation 
in  that  they  are  both  based  on  linearized  ship  wave  theory  and  the  thin 
ship  approximation  in  the  framework  of  the  potential  flow  of  an  incompressibl 
invlscid  fluid.  In  general,  thin-ship  theory  can  be  applied  to  any 
conventional  catamaran  with  arbitrary  camber.  However,  as  this  theory 
has  been  specialized  to  SWATH  configurations,  the  ship  geometry  is 


restricted  to  ships  with  "no-cross  flow"  cambers.  In  practice,  this 
means  SWATHS  with  very  small  strut  thicknesses  compared  to  the  transverse 
strut  separation. 

In  the  theory  of  ship  resistance  extrapolation,  the  resistance  of  a 
ship  is  typically  decomposed  into  residuary  and  frictional  resistance 
components.  The  frictional  resistance  component  is  generally  taken  as 
equal  to  the  resistance  of  a  flat  plate  with  area  of  the  ship,  and  with 
length  equal  to  the  length  of  the  ship.  The  three  dimensional  viscous 
effects  are  then  included  in  the  residuary  resistance  along  with  other 
resistance  components  such  as  wave  making  resistance.  For  a  SWATH  ship, 
the  calculation  of  the  frictional  resistance  is  not  as  simple  as  stated 
above,  however,  the  same  principle  applies  if  one  calculated  the 
frictional  resistance  of  the  bodies  and  struts  separately. 

The  wave  making  resistance  component  of  the  residuary  resistance  can 
be  computed  using  inviscid  fluid  flow  theory.  However,  the  other 
components  of  the  residuary  resistance  are  not  easily  calculated 
theoretically.  Because  of  this,  the  remainder  of  the  residuary  resistance 
is  computed  empirically.  The  difference  between  the  experimentally 
determined  residuary  resistance  and  the  theoretically  determined  wave 
resistance  is  calculated  for  those  cases  where  adequate  experimental  data 
exists.  These  differences  which  are  denoted  as  "form  drag",  are  then 
plotted  as  a  function  of  strut  speed-length  ratio.  Under  the  assumption 
that  the  form  drag  remains  unchanged  for  reasonable  changes  of  SWATH 
proportions,  a  curve  can  be  faired  through  the  individual  form  drag  curves, 
and  applied  to  new  hull  forms. 
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For  the  wave  resistance  predictions,  the  geometries  of  the  body  and 
8trut(s)  are  represented  by  modified  Chebychev  series.  The  methods  used 
to  determine  these  series  are  the  principal  differences  between  the 
programs  CLOSEFIT  and  SYNTHESIS.  The  program  CLOSEFIT  uses  an  integration 
technique,  similar  to  that  used  to  determine  the  coefficients  of  a  Fourier 
series,  to  determine  the  coefficients  of  the  Chebychev  series  based  on 
the  actual  hull  offsets.  The  program  SNYTHESIS  uses  statical  moments  of 
the  strut  waterplane  and  body  to  define  the  coefficients  of  the  Chebychev 
series.  Typically,  the  CLOSEFIT  program  uses  a  series  of  10  to  20  pairs 
of  terms  to  represent  a  body  or  strut,  where  each  pair  of  terms  contains 
an  even  and  an  odd  Chebychev  function,  in  the  same  fashion  that  a  Fourier 
series  consists  of  pairs  of  sine  and  cosine  functions.  The  program 
SYNTHESIS  uses  a  series  consisting  of  only  three  pairs  of  terms  to 
represent  the  body  and  strut.  This  reduction  in  the  number  of  terms 
results  in  savings  of  ten  to  forty  times  in  the  computation  of  the  wave 
resistance  in  SYNTHESIS  as  compared  to  CLOSEFIT. 

In  addition  to  the  effective  power  computations,  the  program  SYNTHESIS 
also  computes  the  delivered  power  which  would  be  required  to  propel  the 
SWATH  ship.  The  propulsion  calculations  are  performed  by  a  subprogram 
SHPCMP  and  its  attendant  subprograms.  SHPCMP  determines  the  wake  and 
thrust  deduction  for  the  design  based  on  empirical  relations  derived  from 
earlier  SWATH  propulsion  experiments.  The  program  then  calculated  the 
optimum  four  bladed  Troost  series  propeller  at  the  design  speed,  along 
with  the  off  design  performance  of  the  propeller.  This  data  along  with 
the  effective  power  computations  is  used  to  calculate  the  delivered  power. 


OVERVIEW  OF  CLOSEFIT 


The  CLOSEFIT  SWATH  program  is  composed  of  several  modules,  each  of 
which  consists  of  a  number  of  subroutines  and  performs  a  specific  func¬ 
tion  in  the  SWATH  ship  resistance  calculation.  Figure  1  presents  a  tree 
diagram  which  reflects  the  relationship  among  all  subroutines. 

Subroutine  CHEB  evaluates  the  coefficients  of  the  Chebychev  series, 
which  represents  the  strut  half  thickness  function  and  the  body 
sectional  area  function  in  wave  resistance  equations.  This  is 
accomplished  by  calling  READ,  SCALE,  COMPUT,  CHECK  and  SURFACE.  The 
offset  tables  are  read  in  by  Subroutine  READ,  then  normalized  by 
Subroutine  SCALE.  The  driver  Subroutine  COMPUT  calls  Subroutine  CHEV, 
which  in  turn  calls  SPLINE  and  SUMSPL  to  compute  the  coefficients  of 
the  Chebychev  series.  The  integrals  used  to  evaluate  the  strut 
Chebychev  coefficients  are  as  follows: 


A 


sm 


=  2  ,  7T/2 
T7  -TT/2 


B 

sm 


=  2.  r  TT/2 
IT  -TT/2 


t(sin  0)  cos(2m-l)6d0 
t  (sin  0)  sin  2m0d0 


where  t(sin  9)  is  a  function  adopted  to  describe  the  offsets  of  the 
strut  or  the  body.  The  Cubic  Spline  Function  Technique  is  used  to  fit  an 
analytic  curve  to  the  discrete  offset  points.  Subroutine  SPLINE  solves 


7 


SWATH 

i 

i 


I — CHE  8 


READ 

■SCALE 

COMPUT- 

CHECK— 

-SURFACE 


CHEV - 

EVAL 

-WTSURFS 

-WTSURFB 


| - SPLINE 

' - SUMSPL 

- EVAL 

- EVAL 


— ECHO 


— PCHEB 


RWAVE 


PLOTI 

-PLOT2 

PLOT3 

-PLOT4- 

-RINIT 

-SIMPSN 

-RINTEG 

-BESSJ 

REFLKT 


QPLOTZ5 


BESSJ 


— PRINT 


> — RWOUT- 


SUM 

FORMDR - YINTP 

CFITTC 


Figure  3  —  Tree  Diagram  of  CLOSEFIT  SWATH  Program 
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the  set  of  simultaneous  linear  algebraic  equations,  which  results  from 
the  spline  fitting  procedure  and  obtains  the  second  derivatives  of  the 
analytical  functions  t(sin  6)  at  each  offset  point.  Next,  Subroutine 
SUMSPL  uses  these  second  derivatives  in  a  set  of  integration  formulas 
based  on  the  Filon  Quadrature  Technique  to  evaluate  the  above  trigono¬ 
metric  integrals.  Subroutine  CHECK  evaluates  the  offsets  approximated 
from  the  Chebychev  series  in  Subroutine  COMPUT  and  prints  out  the 
predicted  offsets  along  with  the  true  offsets.  This  permits  a  check  on 
the  accuracy  of  the  Chebychev  series  approximations. 

Subroutine  SURFACE  computes  the  total  SWATH  ship  wetted  surface  by 
calling  WTSURFS  and  WTSURFB.  Subroutine  ECHO  echoes  the  input  data  of 
the  SWATH  program  on  the  printer  as  a  check  of  its  accuracy.  Subroutine 
PCHEB  plots  the  strut  half  thickness  curve  and  body  sectional  area 
curve  on  the  printer  in  order  to  validate  the  Chebychev  series 
approximation. 

A  large  amount  of  the  execution  time  of  this  program  is  spent  in 
the  evaluation  of  the  auxiliary**  functions ,  T  and  W,  given  by: 
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Due  to  the  highly  oscillatory  nature  of  the  integrands  of  these 
auxiliary  functions,  the  integration  range  is  divided  into  three  separate 
intervals.  For  each  interval,  a  particular  numerical  integration  method 
is  used  to  improve  the  accuracy  of  the  calculation.  To  avoid  the 
singularity  in  the  denominator  of  the  integrand,  the  integration  variable 
(a)  is  replaced  by  Z  =*Ja  -  Yos  in  the  interval  from  yos  to  yQS  +  1. 

The  numerical  method  used  in  this  interval  is  Simpson's  Rule,  which 
uses  30  points.  The  central  region  of  the  integration  process  goes 
from  a  =  yQS  ►  1  to  a  =  a^x.  The  choice  of  ctjnax  is  based  on  the  require¬ 
ment  for  accuracy  in  the  numerical  integration  method  used  in  this  region, 
where  a  step  size  is  computed  at  a  given  a,  and  a  three-point  Simpson's 
Rule  is  used  until  the  upper  integration  limit  of  the  region,  «max.  is 
reached . 

An  analytical  integration  formula  was  developed  to  determine  the 

effects  of  the  final  region,  from  a  =  amax  to  a  =  a£smax  *  wflere  a£  is  an 

~“”'u  smax 

empirical  constant  which  defines  the  upper  limit  of  the  tail  region. 

Beyond  this  point,  the  integration  area  is  assumed  to  be  negligible. 

Subroutine  PRINT  will  print  the  table  of  the  T  and  W  arrays, 
provided  the  flag  JTEST  is  set  other  than  zero. 

Subroutine  RWOUT  calls  SUM  to  evaluate  the  components  of  the  wave 
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resistance,  based  on  Lin  and  Day's  equations^.  Form  drag  is  evaluated 
by  Subroutine  FORMDR,  which  is  based  on  experimental  data.  Frictional 
drag  is  estimated  by  calling  Subroutine  CFITTC,  which  is  baaed  on  the  C 

F 

2 

curve  adopted  by  the  International  Towing  Tank  Conference  (ITTC)  . 

For  each  test  speed,  a  table  is  printed  showing  all  the  components  of 
total  ship  resistance  and  the  effective  power  required  to  tow  the  SWATH 
ship  at  a  constant  speed. 


References  are  listed  on  page  112. 
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OVERVIEW  OF  SYNTHESIS 


The  procedure  of  the  SYNTHESIS  SWATH  program  parallels  the  CLOSEFIT 
program  in  calculating  total  ship  resistance^  except  in  its  approach  to 
obtaining  the  coefficients  of  the  Chebychev  series.  In  SYNTHESIS,  the 
waterplane  area  coefficients,  waterplane  moment,  waterplane  inertia, 
body  prismatic  and  body  moment  are  used  to  evaluate  the  Chebychev 
coefficients.  The  only  other  major  difference  between  the  two 
programs  is  that  SYNTHESIS  includes  a  propeller  design  and  optimization 
program. 

Figure  2  presents  the  tree  diagram  of  all  subroutines  called  by  the 
SYNTHESIS  program.  The  module  i*HEB  is  much  simpler  in  SYNTHESIS  than 
its  CLOSEFIT  counterpart  due  to  its  approach  in  evaluating  the  Chebychev 
coefficients  through  geometric  identities  as  opposed  to  numerical  integrations. 

Modules  RWAVE  and  ROUT  are  identical  for  both  programs,  and  in  both 
require  a  great  deal  of  execution  time. 

Module  FINDRG  evaluates  the  fin  drag  of  a  SWATH  ship.  It  is 
directly  excerpted  from  coding  documented  in  reference  (3). 

Module  SHPCMP  can  determine  propeller  design  and  off-design 
performance  if  the  propeller  diameter  is  specified.  The  subroutines 
used  in  this  module  are  excerpted  from  coding  documented  in  reference  (4) . 
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COMMON/AUX/ 

PURPOSE:  AUX  stores  the  auxiliary  function  table. 


NAME 

TYPE 

LENGTH 

DEFINITION 

TS 

R 

(10,10) 

Auxiliary  function  T 

for 

strut 

WS 

R 

(10,10) 

Auxiliary  function  W 

for 

strut 

TB 

R 

(10,10) 

Auxiliary  function  T 

for 

body 

WB 

R 

(10,10) 

Auxiliary  function  W 

for 

body 

TSB 

R 

(10,10) 

Auxiliary  function  T 
body 

for 

strut 

to 

WSB 

R 

(10,10) 

Auxiliary  function  W 
body 

for 

strut 

to 

TS12 

R 

(10,10) 

Auxiliary  function  T 
to  strut  2 

for 

strut 

1 

WS12 

R 

(8,10,10) 

Auxiliary  function  W 
to  strut  2 

for 

strut 

i 

TB12 

R 

(8,10,10) 

Auxiliary  function  T 
body  2 

for 

body  1  to 

WB12 

R 

(8,10,10) 

Auxiliary  function  W 
body  2 

for 

body 

1  to 

TSB  12 

R 

(8,10,10) 

Auxiliary  function  T 
body  2 

for 

strut 

1  to 

WSB  12 

R 

(8,10,10) 

Auxiliary  function  W 
body  2 

for 

strut 

1  to 

TSBP 

R 

(10,10) 

Auxiliary  function  T 

WSBP 

R 

(10,10) 

Auxiliary  function  W 

TSB12P 

R 

(10,10) 

Auxiliary  function  T 

WSB12P 

R 

(10,10) 

Auxiliary  function  W 

COMMON/NAME/ 


PURPOSE: 

NAME 

INCHAR 


NAME 


NAME  stores  a  flag  indicating  whether  the  offset  data  are 
for  the  strut  or  the  body  and  providing  the  alphanumeric 
information  for  the  strut  and  the  body. 

TYPE  LENGTH  DEFINITION 

I  8  INCHAR  (1)  =  STRUT.  Strut  offset 

data  will  follow. 

INCHAR  (1)  =  BODY.  Body  offset  data 
will  follow. 

INCHAR  (2)  to  INCHAR  (8) .  Hold  the 
alphanumeric  information  for  the  strut 
or  the  body. 

I  7  Temporary  storage  location  for  alpha¬ 

numeric  strut  or  body  information. 
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COMMON/NORM/ 

PURPOSE:  NORM  stores  the  scaling  factors  for  normalizing  X  and  Y. 


NAME 

TYPE 

LENGTH 

DEFINITION 

YMAX 

R 

Maximum  value  of  Y 

XDIVER 

R 

Divisor  for  normalizing  X 

YDIVER 

R 

Divisor  for  normalizing  Y 

XADDER 

R 

Number  to  be  added  to  normalize  X 

) 


COMMON/OMEGA/ 

PURPOSE:  OMEGA  stores  the  variables  used  in  evaluating  the  auxiliary  function. 


NAME 

TYPE 

LENGTH 

DEFINITION 

NMAX 

I 

Maximum  order  of  Chebychev  Series 

HSOLS 

R 

Ratio  of  draft  of  strut  to  length 

of  strut  (h  /L  ) 
s  s 

HBOLB 

R 

Ratio  of  depth  of  submergence  to 
length  of  body  (h^/L^) 

GAMAOS 

R 

vos  -  g  yv2 

GAMAOB 

R 

Yob  •  s  yv2 

NOSEPS 

I 

Number  of  different  hull  separation 
distances 

SEPDIS 

R 

(8) 

Hull  separation  distance  b  (ft) 

SEP 

R 

(8) 

2b/ (Y  -L  ) 

'os  s 

C STRUT 

R 

(8) 

Distance  of  CL  of  strut  to  CL  of  body 

CSTRT2 

R 

(8) 

Distance  of  CL  of  2nd  strut  (if  present) 
from  CL  of  body 

GOSQ 

R 

2 

^os 

h 

PHIS 

R 

2(1T)/Yos 

s 

PHIB 

R 

2<^ob 

RATIOL 

R 

Yob/YoS 

COMMON/OUT/ 


PURPOSE: 

OUT  defines  and 
parameters. 

stores  physical  constants  and  geometric 

NAME 

TYPE 

LENGTH 

DEFINITION  SET 

VALUES 

HS 

R 

Draft  of  strut  h  (ft) 
s 

HB 

R 

Depth  of  submergence  from  surface  to 
body  centerline  h^  (ft) 

XLS 

R 

Length  of  strut  Lg  (ft) 

XLB 

R 

Length  of  body  (ft) 

TSMAX 

R 

Maximum  thickness  of  strut  t  (ft) 

max 

AX 

R 

2 

Maximum  cross-sectional  area  of  body  Ax  (ft  ) 

PI 

R 

3.1415926535897 

G 

R 

2 

Acceleration  due  to  gravity  g  (ft/sec  )  32.155 

RHO 

R 

2  4 

Density  of  water  p  (lb  x  sec  /ft  )  1.9367 

GNU 

R 

2  5 

kinematic  viscosity  v  (ft  /sec)  1.2970  x  10 

WETS 

R 

2 

Wetted  surface  area  of  strut  S  (ft  ) 

s 

WETB 

R 

2 

Wetted  surface  area  of  body  S,  (ft  ) 

D 

WTSURF 

R 

2 

Total  wetted  surface  area  St  (ft  ) 

VMFPS 

R 

Test  velocity  V  (fps) 

DELCF 

R 

Correction  allowance  AC  0.0005 

F 

TITLE 

I 

(8) 

Array  of  characters  containing  the  title 
of  the  test 
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COMMON/PLOT/ 

PURPOSE:  PLOT  stores  the  constants  and  characters  for  the  plotting  routine 


NAME 

TYPE 

LENGTH 

DEFINITION 

SET 

VALUES 

NFIRST 

I 

Position  in  the  arrays  of  the 
ordered  pair  to  be  plotted 

first 

1 

NLA  ST 

I 

Position  in  the  arrays  of  the 
ordered  pair  to  be  plotted 

last 

101 

NPOINT 

I 

NPOINT  equals  "1"  each  point  from 

1 

N FIRST  to  NLA ST  is  to  be  plotted,  "2" 
if  every  other  point  is  to  be  plotted, 
etc. 


NMAX 

R 

Value  of  abscissa  at  right -most  grid 
line 

1.0 

XMIN 

R 

Value  of  abscissa  at  left-most  grid 
line 

-1.0 

NSCLI 

L 

Logical  value  (Should  be  false  if  PLOTI 
has  not  been  called  and  standard  grid  is 
desired) 

NCHAR 

I 

Number  of  valid  characters  in  label 

N SCALE 

I 

(4) 

Printing  scale  factor  of  ordinate 

PCHAR 

I 

(2) 

Plotting  characters 

COMMON/PSI/ 

PURPOSE:  PSI  stores  the  parameters  for  the  Integration  procedure 


NAME 

TYPE 

LENGTH 

DEFINITION 

NPTSZ 

I 

Number  of  integration  steps  from 

Y  to  Y+l  (odd  number) 
os  os 

PTSAF 

R 

Scaling  factor  of  step  size  in 
integrating  from  to  a£gmax 

EXPN 

R 

Empirical  constant  for  integra¬ 
tion  to  stop 

NALMAX 

I 

Maximum  number  of  integration  steps 

from  y+l  to  a 

os  m&x 

NAL 

I 

Counter  of  integration  steps 

TAIL 

R 

Integration  made  from  a  to  a£s 

ALFA 

R 

Integrating  variable  (a) 

ALMAX 

R 

Upper  limit  of  a  for  integration 

N STEPS 

I 

Number  of  integration  steps  from 

a  to  a„ 

max  Slsmax 

COMMON/ XRPLOTF/ 

PURPOSE: 

XRPLOTF 

stores  the  values  of 

variables  for  the  plotting  routine 

NAME 

TYPE 

LENGTH 

DEFINITION 

XL 

R 

Value  of 

abscissa  at  left-most  grid  line 

XH 

R 

Value  of 

abscissa  at  right -most  grid  line 

YL 

R 

Value  of 

ordinate  at  bottom  grid  line 

YH 

R 

Value  of 

ordinate  at  top  grid  line 

XMOV 

R 

Abscissa 

index  increment  numberofor  array 

YMOV 

R 

Ordinate 

index  increment  number  for  array 

GRAF 

GRAF 
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COMMON / XRPLOTG/ 


PURPOSE: 

NAME 

GRAF 


XEPLOTG  stores  the  values  and  characters  of  variables  for  the 


plotting  routein 

TYPE 

LENGTH 

DEFINITION 

I 

(11,204) 

Array  containing  the  image  to  be  plotted 

*  r 


COMMON /XRPLOTQ/ 

PURPOSE:  XRPLOTQ  stores  the  constants  and  characters  for  the  plotting  routine 

SET 

NAME  TYPE  LENGTH  DEFINITION  VALUES 

II  I  Ordinate  scale  factor 

JJ  I  Ordinate  scale  factor 

KK  I  Ordinate  scale  factor 

LL  I  Ordinate  scale  factor 

NHL  I  Number  of  horizontal  grid  lines  6 

NSBH  I  Number  of  spaces  between  adjacent  10 

horizontal  grid  lines 

NVL  I  Number  of  vertical  grid  lines  11 

NSBV  I  Number  of  spaces  between  adjacent  10 

vertical  grid  lines 


NCHAR 


Plotting  character  of  horizontal  grid  1H. 
lines 


VCHAR 


Plotting  character  of  vertical  grid  1H. 
lines 

Number  of  vertical  spaces 


Number  of  horizontal  spaces 

Logical  variable,  *  .  TRUE  .  when  the 
maximum  and  minimum  values  of  the  oridinate 
are  determined 


Logical  variable,  m  .  TRUE  .  when  the 
maximum  and  minimum  values  of  the  abscissa 
are  determined 
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SUBROUTINE  DOCUMENTATION  FOR  CLOSEFIT 


NAME: 


PROGRAM  SWATH 


PURPOSE: 

CALLING  SERVICE: 

ARGUMENTS: 

COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
COMMENTS: 


Program  SWATH  is  the  main  program  for  the 
CLOSEFIT  resistance  prediction  of  Small-Water- 
p lane- Area  Twin-Hull  Ships 

Program  SWATH  (INPUT,  OUTPUT,  TAPE5  -  INPUT 
TAPE6  -  OUTPUT,  TAPE 8) 

NONE 

OUT,  AUX,  OMEGA,  PSI,  COEFS 

CHEB,  ECHO,  PCHEB,  PRINT,  RWAVE,  RWOUT 
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NAME: 


SUBROUTINE  BESSJ 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


SUBROUTINE  BESSJ  evaluates  the  Bessel  function 
from  order  0  to  order  N. 

CALL  BESSJ  (X,  N,  VJ) 

X  Argument  of  the  Bessel  function 
N  Maximum  order  of  the  Bessel  function 
VJ  Array  holding  (N+l)  values  of  the  Bessel 
function  of  order  zero  up  to  N,  where 

VJ (0)  •  J0(X) 


VJ(N)  =  JN(X) 


NONE  \ 

NONE 

RINTEG,  RWAVE 


NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


FUNCTION  CFITTC 

Function  CFITTC  determines  the  ITTC 
frictional  resistance  coefficient 

CF  =  CFITTC  CRN) 

RN  Reynolds  number  at  test  condition 

NONE 

NONE 

RWOUT 

r  -  0.075 

Cf  ~  - - 

flo^0(Rn)  ~2)2 
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NAME: 


SUBROUTINE  CHEB 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  CHEB  computes  from  offset  tables 
the  coefficients  of  the  Chebychev  Series, 
which  is  an  approximation  of  the  outline  of 
the  waterplane  area  or  body  cross-sectional 
area  curve. 

CALL  CHEB  (AMC,  BMC,  MMAX) 

AMC  Coefficients  of  the  Chebychev 

Cosine  Series 

BMC  Coefficients  of  the  Chebychev 

Sine  Series 

MMAX  Maximum  order  of  the  Chebychev 

Series 

NAVME,  OFFSET 

READ,  SCALE,  COMPUT,  CHECK,  SURFACE 
SWATH 
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NAME : 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 


SUBROUTINE  CHEB 

Subroutine  CHEB  computes  from  offset  tables 
the  coefficients  of  the  Chebychev  Series, 
which  is  an  approximation  of  the  outline  of 
the  waterplane  area  or  body  cross-sectional 
area  curve. 

CALL  CHEB  (AMC,  BMC,  M4AX) 

AMC  Coefficients  of  the  Chebychev 

Cosine  Series 

BMC  Coefficients  of  the  Chebychev 

Sine  Series 

MMAX  Maximum  order  of  the  Chebychev 

Series 

NAME,  OFFSET 

READ,  SCALE,  COMPUT,  CHECK,  SURFACE 
SWrtTH 


COMMENTS : 


NAME: 


SUBROUTINE  CHECK 


PURPOSE:  Subroutine  CHECK  evaluates  the  Chebyehev  Series 


from  the  Chebyehev  coefficients  obtained  from 
Subroutine  COMPUT  and  prints  out  the 
evaluated  offset  table  for  manual  checking. 

CALLING  SEQUENCE: 

CALL  CHECK  (AMC,  BMC,  MMAX,  ISBODY) 

ARGUMENTS : 

AMC 

Coefficients  of  the  Chebyehev 

Cosine  Series 

BMC 

Coefficients  of  the  Chebyehev 

Sine  Series 

MMAX 

Maximum  order  of  the  Chebyehev 
polynomial 

ISBODY 

Logical  flag,  indicating  "Is  this 
for  body?" 

i  .  TRUE  .  for  body 

ISBODY 

= 

1  .  FALSE  .  for  strut 

COMMON  BLOCKS: 

OFFSET, 

NAME 

SUBROUTINE  CALLED: 

EVAL 

CALLED  BY: 

CHEB 

COMMENTS : 


NAME: 


SUBROUTINE  CHEV 


PURPOSE:  Subroutine  CHEV  uses  the  offsets  to  compute 

the  Chebychev  coefficients  via  spline 
approximation  and  numerical  integration. 

CALLING  SEQUENCE:  AMC  Coefficients  of  the  Chebychev 

Cosine  Series 

BMC  Coefficients  of  the  Chebychev 

Sine  Series 

MMAX  Maximum  order  of  the  Chebychev 

polynomial 

ISBODY  Logical  flag,  indicating  "Is  this 
for  body?" 
f .  TRUE  .  for  body 
ISBODY  = < 

1^.  FALSE  .  for  strut 
COMMON  BLOCKS:  OFFSET,  NAME 

SUBROUTINES  CALLED:  SPLINE,  FUNCTION  SUMPSPL 

CALLED  BY:  COMPUT 

COMMENTS : 


NAME: 


SUBROUTINE  COMPUT 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


Subroutine  COMPUT  initiates  the  computation 

of  the  Chebychev  coefficients. 

CALL  COMPUT  (AMC,  BMC,  MMAX,  ISBODY) 

AMC  Coefficients  of  the  Chebychev 

Cosine  Series 

BMC  Coefficients  of  the  Chebychev 

Sine  Series 

MMAX  Maximum  order  of  the  Chebychev 

polynomials 

ISBODY  Logical  flag,  indicating  "Is  this 

for  body?" 

f .  TRUE  .  for  body 

ISBODY  H 

I  .  FALSE  .  for  strut 

v_ 

OFFSET,  NAME,  NORM 

CHEV 

CHEB 


NAME: 


SUBROUTINE  ECHO 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  ECHO  prints  a  summary  of  input  data  and 
the  Chebychev  coefficients 

CALL  ECHO  (NLOC) 

NLOC  Number  of  strut  locations  to  be 

tested 

AUX,  OMEGA,  OUT 

NONE 

SWATH 
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NAME: 


SUBROUTINE  EVAL 


PURPOSE:  Subroutine  EVAL  evaluates  the  Chebychev  Series 

at  a  given  station. 


CALLING  SEQUENCE: 

CALL  EVAL 

(XX,  MAX,  FFSET,  AMC,  BMC) 

ARGUMENTS : 

XX 

Value  of  X  station  on  a  range  of 
-1  to  +1 

MMAX 

Maximum  order  of  the  Chebychev 
polynomial 

FFSET 

Offset  approximated  by  the  Chebychev 
Series 

AMC 

Coefficients  of  the  Chebychev  Cosine 
Series 

BMC 

Coefficients  of  the  Chebychev  Sine 
Series 

COMMON  BLOCKS: 

NONE 

SUBROUTINE  CALLED: 

NONE 

CALLED  BY: 

CHECK 

COMMENTS : 


i 


( 

i 

I' 


1 


NAME: 


FUNCTION  FORMDR 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 


Function  FORMDR  evaluates  the  form 

drag  coefficient,  based  on  experimental  data 

FOR  =  FORMDR  (VL) 

VL  Ratio  of  ship  speed  to  length  of 

strut  V//CJT 

NONE 

FUNCTION  YINTP 
ROUT 

_3 

Form  drag  Is  assumed  to  be  0.5  x  10 

for  less  than  1.325  and  greater  than  2.5 


! 


COMMENTS: 


NAME: 


SUBROUTINE  PCHEB 


PURPOSE:  Subroutine  PCHEB  plots  by  line  printer  the 

body  sectional  area  curve  and  waterplane 
outline  curve  from  the  given  Chebychev 
coefficients 

CALLING  SEQUENCE:  CALL  PCHEB  (AS,  BS,  AB,  BB,  NN,  TITLE) 

ARGUMENTS:  AS  Coefficients  of  Chebychev  Sine 

Series  for  strut 

BS  Coefficients  of  Chebychev  Cosine 

Series  for  strut 

AB  Coefficients  of  Chebychev  Sine 

Series  for  body 

BB  Coefficient  of  Chebychev  Cosine 

Series  for  body 

NN  Maximum  order  of  Chebychev  Series 

TITLE  Array  containing  the  alphanumeric 

characters  of  the  title  of  the  experi¬ 
ment 

COMMON  BLOCKS:  PLOT,  XRPLOTQ 

SUBROUTINES  CALLED:  PL0T1,  PLOT 2,  PL0T3 ,  PLOT 4 

CALLED  BY:  CHEB 

COMMENTS: 
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NAME: 


SUBROUTINE  PLOT1 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCK: 
SUBROUTINE  CALLED: 
CALLED  BY 
COMMENTS: 


Subroutine  PLOT1  sets  up  spacing  and 
determines  the  values  of  the  axes. 


CALL  PLOT  1  (NSCALE ,  A,  B,  C,  D,  E,  F) 


NSCALE 

NSCALE  (1) 

NSCALE  (2) 

NSCALE  (3) 

NSCALE  (4) 


E 

F 

XRPLOTQ 


Integer  array  defined  as 
follows: 

I,  if  printed  values  of  the 
ordinate  are  10  **  I  times 
the  actual  value 

J,  if  printed  values  of  the 
ordinate  are  10  **  J  times 
the  actual  value 

K,  if  printed  values  of  the 
abscissa  are  10  **  K  times 
the  actual  values 

L,  if  printed  values  of  the 
abscissa  are  10  **  L  times 
the  actual  value 

Integer  number  of 
horizontal  grid  lines 

Integer  number  of  spaces 
beyond  each  horizontal  grid 
line  to  the  next  grid  line 

Integer  number  of  vertical 
grid  lines 

Integer  number  of  spaces 
beyond  each  vertical  grid 
line  to  the  next  grid  line 

Horizontal  grid  character 

Vertical  grid  character 


NONE 


PCHEB 
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NAME: 


SUBROUTINE  PLOT2 


PURPOSE: 


CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


Subroutine  PLOT2  examines  the  minimum  and 
maximum  values  of  the  abscissa  and  the 
ordinate  and  establishes  an  internal 
formula  for  computing  location  in  the  image 
region  corresponding  to  the  point  to  be 
plotted. 

CALL  PLOT 2  (XMAX,  XMIN,  YMAX ,  YMIN,  NSCLI) 

XMAX  Value  of  abscissa  at  right-most 

grid  line 

XMIN  Value  of  abscissa  at  left-most 

grid  line 

YMAX  Value  of  ordinate  at  top  grid 

line 

YMIN  Value  of  ordinate  at  bottom  grid 

line 

NSCLI  Logical  flat  (should  be  .FALSE., 

if  PL0T1  has  not  been  called  and 
standard  grid  is  desired) 

XRPLOTF ,  SRPLOTQ,  XRPLOTG 

NONE 

PCHEB 
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NAME: 


SUBROUTINE  PLOT3 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  PLOT3  assigns  an  alpha-character 
to  each  point  to  be  plotted. 

CALL  PLOT3  (PCHAR,  X,  Y,  SDATA,  FDATA ,  DDATA) 

PCHAR  Plotting  character 

X  Array  containing  the  X  coordinates 

to  be  plotted 

Y  Array  containing  the  Y  coordinates 

to  be  plotted 

SDATA  Integer  position  in  the  arrays  of 

the  first  ordered  pair  to  be  plotted 
Cl  if  each  point  from  SDATA  to  DDATA 
j  is  to  be  plotted 

FDATA  =<,  2  if  every  other  point  is  to  be 
plotted 

3  if  every  third  point  is  to  be 
^  plotted 

DDATA  Integer  position  in  the  array  of  the 
last  ordered  pair  to  be  plotted 

XRPLOTF,  XRPLOTG 

NONE 

PCHEB 
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NAME: 


SUBROUTINE  PL0T4 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


Subroutine  PLOT4  prints  the  image  of  the 
completed  graph  on  the  printer,  including 
the  values  of  the  abscissa  and  the  ordinate 
at  the  grid  lines  outside  the  bottom  and 
left  edge  of  the  graph. 

CALL  PLOT 4  (MCHAR,  NCHAR) 

MCHAR  Single  dimension  array  containing 

alpha-characters  to  be  plotted  at 
the  left  of  the  graph 
NCHAR  Number  of  valid  characters  in 

MCHAR 

XRPLOTF,  XRPLOTG,  XRPLOTQ 

QPL0TZ5 

PCHEB 


NAME: 


SUBROUTINE  PRINT 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  PRINT  prints  the  auxiliary 
function  tables  of  T  and  W. 


CALL  PRINT  (I,  NLOC2) 


I 

Index  indicating  strut  location 

NLOC2 

Flag  indicating  the  presence  of 
a  second  strut 

| 

0  if  single  strut 

NLOC2  = 

1  if  tandem  struts 

OUT,  AUX, 

OMEGA 

NONE 


SWATH 


NAME: 


SUBROUTINE  QPLOTZ5 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
COMMENTS : 


Subroutine  QPLOTZ5  calculates  the  scaling 
information  needed  to  generate  the  format 
to  label  the  left-hand  side  of  the  program. 

CALL  QPLOTZ5  (PDQ) 

PDQ  Scaling  factor  for  ordinate  plot 

XRPLOTF ,  XRPLOTQ 

NONE 
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NAME: 


SUBROUTINE  READ 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


Subroutine  READ  reads  the  offset  data  of 
the  strut  or  body  and  checks  for  error. 

CALL  READ 

NONE 

OFFSET,  NAME,  NORM 

NONE 

CHEB 
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NAME: 


SUBROUTINE  REFLKT 


V 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS : 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  REFLKT  reflects  all  the 
symmetrical  matrices  of  T  and  W. 

CALL  REFLKT 

NONE 

AUX,  OMEGA 

NONE 

RWAVE 


1 
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NAME: 

PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


SUBROUTINE  RINIT 

Subroutine  RINIT  initializes  the  T  and  W 
arrays  to  zero. 

CALL  RINIT 

NONE 

AUX,  OMEGA 

NONE 

RWAVE 
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NAME: 


SUBROUTINE  RINTEG 


PURPOSE : 


Subroutine  RINTEG  evaluates  the  integrand 
for  T  and  W  functions. 


CALLING  SEQUENCE: 

ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 


CALL  RINTEG  (ALFA,  B,  D,  NLOC2,  WTINT, 
SEPCOS,  SQ) 


ALFA 

B 

D 

NLOC2 


NLOC2  = 


WTINT 


SEPCOS 


SQ 


Integrating  variable  a 

cstrut(i)/xls 

CSTRT2(I)/XLS 

Flag  indicating  the  presence  of 
a  second  strut 
'O  if  single  strut 

1  if  tandem  struts 
“Weighting  constant  for  the 
integrand 

Value  of  the  cosine  function  in 
the  integrand 

Value  of  part  of  the  integrand 


(a2 


AUX,  OMEGA 
BESSJ 


CALLED  BY: 


RWAVE 


COMMENTS : 


NAME: 


SUBROUTINE  RWAVE 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS: 


Subroutine  RWAVE  computes  the  auxiliary 
function  of  T  and  W. 


CALL  RWAVE  (B,  D,  NL0C2) 


B 

D 

NLOC2 


NL0C2 


CSTRUT  (I)/XLS 
CSTRT2  (I)/XLS 

Flag  indicating  the  presence 
of  a  second  strut 
-0  if  single  strut 

■1  if  tandem  struts 


OUT,  AUX,  OMEGA,  PS I 

RINIT,  SIMPSN,  RINTEG,  BESSJ,  REFLKT 

SWATH 


Typical  Auxiliary  Functions 


Smw 


CmH2n)  ) 


f 

Tos 


da 


1  s/o1  --1 


J- <“> 
2b  <“> 


where 
D 


E,  -  |  -  e-»<V-s><<»}/iW 


>01“ 


«LS 

<2U*> 
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NAME: 


SUBROUTINE  RWOUT 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  RWOUT  computes  and  prints  drag 
coefficients  and  total  resistance  of  the 
test  model. 

CALL  RWOUT  (I,  NLOC2) 

I  Index  indicating  strut  location 

NLOC2  Flag  indicating  the  presence  of 

a  second  strut 

{0  if  single  strut 

1  if  tandem  struts 
OUT,  AUX,  OMEGA 

SUM,  FUNCTION  FORMDR ,  FUNCTION  CFITTC 
SWATH 
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NAME: 


SUBROUTINE  SCALE 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  SCALE  computes  scaling  factors 
and  scales  the  input  for  length  between 
-1.0  and  +1.0  maximum  beam  and  sectional 
area  to  1.0. 

CALL  SCALE  (ISBODY) 

ISBODY  Logical  flag  indicating,  "Is  this 
for  body?" 
f .  TRUE  .  for  body 
ISBODY  =  <^ 

(_.  FALSE  .  for  strut 
OFFSET,  NAME,  NORM 
NONE 
CHEB 


< 


NAME: 


SUBROUTINE  SIMPoN 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  SIMPSN  sets  up  Simpson's 
multipliers  for  numerical  integration  by 
Simpson's  rule. 

CALL  SIMPSN  (NPTS ,  SIMP) 

NPTS  Number  of  integration  steps 

SIMP  Array  containing  Simpson's 

multipliers 

NONE 

NONE 

RWAVE 
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NAME: 


SUBROUTINE  SPLINE 


PURPOSE: 

CALLING  SEQUENCE: 

ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 


Subroutine  SPLINE  fits  smooth  spline 
segments  through  a  given  set  of  discrete 
data  points. 


CALL  SPLINE  (X,  Y,  D,  N,  BC,  KODE,  KKODE , 
ISBODY) 


X  Value  of  abscissas  of  offset  data 
Y  Value  of  ordinates  of  offset  data 
D  Value  of  second  derivatives  at  a 


N 

BC 

KODE  = 


KKODE  = 
ISBODY 

ISBODY  = 


r 


1 


point 

Total  number  of  offset  data 
points 

Value  of  boundary  conditions 
0  if  user  is  specifying  boundary 
conditions 

1  if  user  is  taking  an  extra¬ 
polation  of  second  derivatives 
as  boundary  conditions 
0  if  there  is  no  printout  from 
spline  routine 

1  if  spline  prints  abscissas, 
ordinates  and  second  derivatives 
Logical  flag  indicating,  "Is  this 
for  body?" 

•  TRUE  •  for  b°dy 


FALSE  •  for  strut 


NONE 


NONE 


CHEV 


COMMENTS: 


1 


NAME: 
PURPOSE : 


CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


SUBROUTINE  SUM 

Subroutine  SUM  computes  sums  of  the  form: 


M  N 

^  ^  {^sm  Asn  Tsmn  +  Bsm  Bsn  Wsmn} 

m=l  n=l 

CALL  SUM  CL,  SUM1S,  SUM1B,  SUM1SB,  SUM12S, 
SUM12B,  SUM12SB) 


L 

Number  of  different  hull 
separation  distances 

SUM1S 

Partial  sum  for  strut  1 

SUM1B 

Partial  sum  for  body  1 

SUM1SB 

Partial  sum  for  interaction 
between  strut  1  and  body  1 

SUM12S 

Partial  sum  for  interaction 
between  strut  1  and  strut  2 

SUM12B 

Partial  sum  for  interaction 
between  body  1  and  body  2 

SUM12SB 

Partial  sum  for  interactions 
between  strut  1  and  body  2  or 
strut  2  and  body  1 

AUX,  OMEGA,  COEFS 


NONE 

RWOUT 


SUM1S 


SUH1B 
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NAME: 


FUNCTION  SUMSPL 


PURPOSE:  Function  SUMSPL  uses  the  spline  coefficients 

and  integrates  the  trigonometric  integrals 
F(X)*SIN(AK*X)  or  F(X)*COS(AK*X) . 

CALLING  SEQUENCE:  SUM  =  SUMSPL  (X,  Y,  S,  N,  AK,  KODE) 


ARGUMENTS:  X 

Y 

S 

N 


AK 

KODE 

COMMON  BLOCKS: 

NONE 

SUBROUTINE  CALLED: 

NONE 

CALLED  BY: 

CHEV 

Values  of  abscissas  of  offset  data 
Values  of  ordinates  of  offset  data 
Values  of  second  derivatives 
Total  number  of  data  points 
(2m- 1)  or  (2m\  m  =  1,2,3  .  .  .  MMAX 
f  +1  if  sine  integration  desired 

t_  -1  if  cosine  integration  desired 


NAME: 

PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS: 


SUBROUTINE  SURFACE 

Subroutine  SURFACE  finds  the  wetted  surface 

area  of  the  strut  or  body  as  well  as  the 

waterplane  area  or  displaced  volume. 

CALL  SURFACE  (AMC,  BMC,  MMAX,  ISBODY) 

AMC  Coefficients  of  the  Chebychev 

Cosine  Series 

BMC  Coefficients  of  the  Chebychev 

Sine  Series 

MMAX  Maximum  order  of  the  Chebychev 

polynomials 

ISBODY  Logical  flag  indicating  "Is  this 

for  body?" 

{•  TRUE  •  for  body 

.  FALSE  .  for  strut 

OUT,  NAME 
WTSURFS ,  WTSURFB 
CHEB 


57 


r 

NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCK: 
SUBROUTINE  CALLED: 
CALLED  BY: 


SUBROUTINE  WTSURFB 

Subroutine  WTSURFB  determines  the  surface 
area  of  a  body  of  revolution  given  by  the 
Chebychev  coefficients. 

CALL  WTSURFB  (AMC,  BMC,  MMAX) 

AMC  Coefficients  of  the  Chebychev 

Cosine  Series 

BMC  Coefficients  of  the  Chebychev 

Sine  Series 

MMAX  Maximum  order  of  the  Chebychev 

Series 

OUT 

EVAL 

SURFACE 


COMMENTS : 


NAME: 
PURPOSE : 


CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCK: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


SUBROUTINE  WTSURFS 

Subroutine  WTSURFS  determines  the  wetted 
surface  area  of  a  strut  whose  thickness 
distribution  is  given  by  the  Chebyehev 
coefficients. 

CALL  WTSURFS  (AMC,  BMC,  MMAX) 

AMC  Coefficients  of  Chebyehev  Cosine 

Series 

BMC  Coefficients  of  Chebyehev  Sine 

Series 

MMAX  Maximum  order  of  Chebyehev 

Series 

OUT 

EVAL 

SURFACE 
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NAME: 


FUNCTION  YINTP 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


Function  YINTP  interpolates  through  a  set  of 
discrete  data. 

FR  =  YINTP  (XA,  X,  Y,  N) 

XA  Point  to  be  interpolated 

X  Array  of  ordinate  data 

Y  Array  of  abscissa  data 

N  Number  of  data  points 

NONE 

NONE 

FORMDR 

Uses  linear  interpolation 
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PROGRAM  DOCUMENTATION  OF  SYNTHESIS 
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COMMON/ AUX/ 

PURPOSE:  AUX  stores  the  auxiliary  wave  resistance  function  table. 


NAME 

TYPE 

LENGTH 

DEFINITION 

TS 

R 

(3,3) 

Auxiliary  function  T 

for 

strut 

WS 

R 

(3,3) 

Auxiliary  function  W 

for 

strut 

TB 

R 

(3,3) 

Auxiliary  function  T 

for 

body 

WB 

R 

(3,3) 

Auxiliary  function  W 

for 

body 

TSB 

R 

(3,3) 

Auxiliary  function  T 
body 

for 

strut  to 

WSB 

R 

(3,3) 

Auxiliary  function  W 
body 

for 

strut  to 

TS12 

R 

(3,3) 

Auxiliary  function  T 
to  strut  2 

for 

strut  1 

WS12 

R 

(3,3) 

Auxiliary  function  W 
to  strut  2 

for 

strut  1 

TB12 

R 

(3,3) 

Auxiliary  function  T 
body  2 

for 

body  1  to 

WB12 

R 

(3,3) 

Auxiliary  function  W 
body  2 

for 

body  1  to 

TSB12 

R 

(3,3) 

Auxiliary  function  T 
body  2 

for 

strut  1  to 

WSB12 

R 

(3,3) 

Auxiliary  function  W 
body  2 

for 

strut  1  to 

TSBP 

R 

(3,3) 

Auxiliary  function  T 

WSBP 

R 

(3,3) 

Auxiliary  function  W 

TSB12P 

R 

(3,3) 

Auxiliary  function  T 

WSB12P 

R 

(3,3) 

Auxiliary  function  W 

> 

*  + 


COMMON/BARCK/IE  / 


PURPOSE: 

NAME 

IE 


IE  = 


BARCK  stores  the  variables  which  define  the  propulsor 
design  condition. 

TYPE  LENGTH  DEFINITION 

I  Error  message  pointer 


r0,if  program  ran  in  normal  fashion 

100, if  limiting  value  of  P/D  has  been  used 

200, if  there  is  an  iteration  malfunction  —  PD(2) 
used  for  P/D 

\  3, if  maximum  BAR,  0.95,  used  --  marine  screw 

4, if  minimum  BAR,  .4,  used  --  marine  screw 

5, if  maximum  tip  speed  used  --  air  screw 

6, if  blade  number  is  outside  program  limit  -- 
V__  closest  number  used 


i 


COWON/CDESC/ 


PURPOSE:  CDESC  stores  the  values  of  constants  which  define  water 

properties . 


NAME 

TYPE 

LENGTH 

DEFINITION 

RHOD 

R 

(3) 

Water  density  array 

ANUD 

R 

(3) 

Kinematic  viscosity  array 

IP 

I 

Pointer  indicating  the  particular 
water  condition  used.  Set  IP  =  1 
for  this  program  to  obtain  sea  water 
condition. 

I 


CQMHON/COEFS/ 

PURPOSE: 

COEFS  stores  the  values  of  the  Chebyehev  coefficients  for 
the  strut  and  body. 

NAME 

TYPE 

LENGTH 

[  DEFINITION 

ASM 

R 

(3) 

Coefficients 

of  Chebychev  Sine  Series  for  strut 

BSM 

R 

(3) 

Coefficients 
for  strut 

of  Chebychev  Cosine  Series 

ABM 

R 

(3) 

Coefficients 

body 

of  Chebychev  Sine  Series  for 

BBM 

R 

(3) 

Coefficients 
for  body 

of  Chebychev  Cosine  Series 

MMAX 


I 


Maximum  order  of  Chebychev  Series 


COMMON/CPROP/ 


PURPOSE: 

CPROP 

used 

stores  the  values  of  the  variables  and  constants 
in  the  propulsor  design  routine. 

NAME 

TYPE 

LENGTH 

DEFINITION 

IPTYP 

I 

(3) 

Indicating  propeller  type 

DIAM 

R 

(3) 

Diameter  of  propeller  (ft) 

VADES 

R 

(3) 

Speed  of  advance  of  propulsor  at  design 
condition  (kts) 

RPMDES 

R 

(3) 

RPM  of  propeller  at  design  point 

TDES 

R 

(3) 

Thrust  of  propulsor  system  (lbs) 

QDES 

R 

(3) 

Torque  of  propulsor  system  (ft-lbs) 

EFFDES 

R 

(3) 

Efficiency  of  propulsor  at  design 
condition 

OMWTD 

R 

(3) 

Thrust  wake  (l-wT) 

CAV 

R 

(3) 

Percent  back  cavitation  allowed  on  water 
screw 

Z 

I 

(3) 

Number  of  blades  of  propeller 

H 

R 

(3) 

Depth  of  water  to  centerline  of 
propulsor  shaft 

BAR 

R 

(3) 

Blade  area  ratio 

PD 

R 

(3) 

P/D  of  propulsor 

J 

R 

(3) 

J  of  propulsor 

I ERROR 

I 

(3) 

Error  message  indicator 

KT 

R 

(3) 

Thrust  coefficient  (K^) 

KQ 

R 

(3) 

Torque  coefficient  (K^) 

COMMON/INPUT/ 

PURPOSE:  INPUT  stores  input  data. 


NAME 

TYPE 

LENGTH 

DEFINITION 

XLS 

R 

Length  of  strut  Lg  (ft) 

HS 

R 

Draft  of  strut  Hg  (ft) 

TSMAX 

R 

Maximum  thickness  of  strut  t  (ft) 

max 

CWP 

R 

Waterplane  area  coefficient  C^,  -  A„/(L-  t 
where  Ay  ■  waterplane  are  of  one  strut  max 

CLCF 

R 

Waterplane  moment  coefficient 

9lcf  ■  v'vv* 

CIYY 

R 

Waterplane  inertia  coefficient, 

ctw  *  1  •  K  >* 

IYY  wp  w  S 

XLB 

R 

Length  of  body  L_  (ft) 

□ 

HB 

R 

Draft  of  body  H  (ft) 
b 

AX 

R 

2 

Maximum  cross-sectional  area  of  body  (ft  ) 

CP 

R 

Body  prismatic  coefficient  Cp  =  ^./(A^  L.  ) 
where  Vj,  “  displaced  volume  of  one  body 

CLCB 

R 

Body  moment  coefficient 

CLCB  "  •  Le)* 

SEPDIS 

R 

Separation  distance  between  centerline  of  the 
two  bodies  (ft) 

CSTRUT 

R 

Distance  from  centerline  of  strut  to 
centerline  of  body  (ft) 

CSTRT2 

R 

Distance  from  centerline  of  second  strut 
to  centerline  of  body  (ft) 

PDIA 

R 

Diameter  of  propeller  D  (ft) 

SPAN 

R 

Span  of  a  demi-fin  (ft) 

CHORD 

R 

Chord  of  a  fin  (ft) 

TFINS 

R 

Maximur.  thickness  of  a  demi-fin  (ft) 

NIOC 

I 

Switch  indicating  the  presence  of  a 
second  strut 

*A11  moments  and  moments  of  inertia  arc  taken  about  the  mid-length  of  the 
respective  strut  or  body. 
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COMMON/OFFDES/ 

PURPOSE:  OFFDES  stores  the  variables  which  define  the  coefficients 

of  KT/KQ  curves. 

NAME  TYPE  LENGTH  DEFINITION 

PJAY  R  (4)  Coefficients  of  polynomial  representing 

KT  or  KQ. 


COM40N/0MEGA/ 


PURPOSE: 

OMEGA  stores  the  values  of  variables  and  constants  for 
evaluation  of  auxiliary  functions  and  ship  resistance 
coefficients . 

NAME 

TYPE  LENGTH 

DEFINITION 

VMFPS 

R 

Speed  of  ship  v  (kts) 

GAMAOS 

R 

2 

Yos  “  **  8  Lg/V  where  U  is  speed  in  fps 

GAMAOB 

R 

Yob  =  g  1^/V2 

GOSQ 

R 

(y  ) 2 

v  'OS'* 

HSOLS 

R 

Ratio  of  draft  to  length  of  strut  H  /L 

j  2) 

HBOLB 

R 

Ratio  of  draft  to  length  of  body  H^/L^ 

WETS 

R 

2 

Wetted  surface  area  of  strut  S  (ft  ) 

s 

WETB 

R 

Wetted  surface  area  of  body  s  (ft2) 

b 

WETFIN 

R 

Wetted  surface  area  of  a  fin  (ft2) 

WTSURF 

R 

Total  wetted  surface  area 

SEP 

R 

2b*W  Lg) 

PHIS 

R 

PHIB 

R 

RATIOL 

R 

CFS 

R 

Frictional  drag  coefficient  of  strut  C_ 

CFB 

R 

s 

Frictional  drag  coefficient  of  body  C_, 

Appendage  drag  (lb) 


APPDRG  R 


COMMON/PCOF  / 

PURPOSE:  PCOF  stores  coefficients  for  KQ/KT  polynomial 

approximation. 

NAME  TYPE  LENGTH  DEFINITION 

CT  R  (4)  Coefficients  of  the  polynomial  which 

is  an  approximation  of  the  KT  curve 

CQ  R  (4)  Coefficients  of  the  polynomial  which 

i<?  an  aoDroximation  of  the  KQ  curve 


COMMON/ PHYSCO/ 

PURPOSE:  PHYSCO  stores  physical  constants. 


NAME 

TYPE 

LENGTH 

DEFINITION 

RHO 

R 

Density  of  water 

GNU 

R 

Kinematic  viscosity  of  water 

G 

R 

Acceleration  due  to  gravity 

PI 

R 

Ratio  of  circumference  to  diameter  of  a 
circle 

DELCF 

R 

Correction  allowance  of  friction  drag 
from  ITTC  line. 

2 


COMMON /PLOT/ 


PURPOSE:  PLOT  stores  the  constants  and  characters  for  the  plotting 
routine . 


NAME 

TYPE 

LENGTH 

DEFINITION 

NFIRST 

I 

Position  in  the  arrays  of  the  first 
ordered  pair  to  be  plotted 

NLAST 

I 

Position  in  the  arrays  of  the  last 
ordered  pair  to  be  plotted 

NPOINT 

I 

NPOINT  equals  "1"  if  each  point  from 
NFIRST  to  NLAST  is  to  be  plotted,  "2" 
if  every  other  point  is  to  be  plotted, 
etc. 

XMAX 

R 

Value  of  abscissa  at  right-most  grid  lin< 

XMIN 

R 

Value  of  abscissa  at  left-most  grid  line 

NSCLI 

L 

Logical  value  (Should  be  false  if  PL0T1 
has  not  been  called  and  standard  grid  is 
desired) 

NCHAR 

I 

Number  of  valid  characters  in  label 

NSCALE 

I 

(4) 

Printing  scale  factor  of  ordinate 

PCHAR 

I 

(2) 

Plotting  characters 

COMMON/PSI/ 

PURPOSE:  PSI  stores  the  values  of  the  variables  and  constants  used 

in  integration  routines. 


LENGTH 


DEFINITION 


NPTSZ 


PTSAF 


NALMAX 


Number  of  integration  steps  from  y  to  y+1 

Scaling  factor  of  step  size  in  integrating 
from  ctjnax  to  ajlsmax 

Empirical  constant  for  integration  to  stop 

Maximum  number  of  integration  steps  from 
Y+1  to 


ALSMAX 


NSTEPS 


Counter  of  integration  steps 

Integration  correction  made  from  to 

aismax 

Integrating  variable  (a) 

Maximum  of  a  for  integration 

Number  of  integration  steps  from  amax  to 
a£smax 


COMMON/ ROOTC/ 


PURPOSE:  ROOTC  stores  the  coefficients  of  a  third  order  polynomial 

which  is  the  approximation  of  the  KT  curve  of  a  propeller. 

NAME  TYPE  LENGTH  DEFINITION 

A  R  4  Array  holding  the  coefficient  of  a  third 

order  polynomial  which  is  the  approximation 
of  the  KT  curve  of  a  propeller 


COMMON / XRP  LOTF / 


PURPOSE: 

XRP LOTF  stores 
routine. 

the  values  of  variables 

for  the  plotting 

NAME 

TYPE  LENGTH 

DEFINITION 

XL 

R 

Value  of  abscissa  at 

left -most  grid  line 

XH 

R 

Value  of  abscissa  at 

right -most  grid  line 

YL 

R 

Value  of  ordinate  at 

bottom  grid  line 

YH 

R 

Value  of  ordinate  at 

top  grid  line 

XMOV 

R 

Abscissa  index  increment  number  for  array 

GRAF 

YMOV 

R 

Ordinate  index  increment  number  for  array 

GRAF 
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COMMON/XRPLOTG/ 

PURPOSE:  XRPLOTG  stores  the  values  and  characters  of  variables  for 

the  plotting  routine. 


NAME 

GRAF 


TYPE  LENGTH  DEFINITION 

I  (11,204)  Array  containing  the  image  to  be  plotted 


1 
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COMMON/XRPLOTQ 


PURPOSE 

:  XRPLOTQ  stores 

routine. 

the  constants  and  characters  for  the  plotting 

NAME 

TYPE 

LENGTH 

DEFINITION 

I 

I 

Ordinate  scale  factor  is  10* 

J 

I 

Number  digits  following  ordinate  decimal 
point 

K 

I 

Abscissa  scale  factor  is  10K 

L 

I 

Number  digits  following  abscissa  decimal 
point 

NHL 

I 

Number  of  horizontal  grid  lines 

NSBH 

I 

Number  of  spaces  between  adjacent 
horizontal  grid  lines 

NVL 

I 

Number  of  vertical  grid  lines 

NSBV 

I 

Number  of  spaces  between  adjacent 
vertical  grid  lines 

HCHAR 

I 

Plotting  character  of  horizontal 
grid  lines 

ISX 

I 

Number  of  horizontal  spaces 

ISY 

I 

Number  of  vertical  spaces 

V 

L 

Logical  variable,  =  .  TRUE  .  when  the 
maximum  and  minimum  values  of  the  ordinate 
are  determined 

H 

L 

Logical  variable,  =  .  TRUE  .  when  the 
maximum  and  minimum  values  of  the  abscissa 
are  determined 
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SUBROUTINE  DOCUMENTATION  FOR  SYNTHESIS 


✓ 
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NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 


COKWON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


SUBROUTINE  BESSJ 

Subroutine  BESSJ  evaluates  the  Bessel  function 
from  order  0  to  order  N. 

CALL  BESSJ  (X,  N,  VJ) 

X  Argument  of  the  Bessel  function 
N  Maximum  order  of  the  Bessel  function 
VJ  Array  holding  (N+l)  values  of  the 
Bessel  function  of  order  zero  up  to 
N,  where 

VJ(0)  =  J0(X) 


VJ (N)  =  JN(X) 

NONE 

NONE 

RINTEG,  RWAVE 


i; 
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NAME: 


SUBROUTINE  BURIL 


PURPOSE:  SUBROUTINE  BURIL  uses  Burrill's  cavitation 

criteria  to  determine  if  the  blade  area 
ratio  (BAR)  used  is  sufficient  and,  if  it 
is  not,  to  estimate  what  BAR  is  necessary. 

CALLING  SEQUENCE:  CALL  BURIL  (T,  D,  H,  VA,  N,  CAVI ,  BAR, 

I ERROR,  PD,  NRET) 


ARGUMENTS : 

T 

D 

H 

VA 

N 

CAVI 

BAR 

I ERROR 

PD 

NRET 

NRET 

CO»WON  BLOCK: 

BARCK 

SUBROUTINE  CALLED: 

NONE 

CALLED  BY: 

PRODES 

Thrust  of  propeller  (lbs) 

Diameter  of  propeller  (ft) 

Depth  of  propeller  shaft  (ft) 

Speed  of  advance  (fps) 

RPM  of  propeller 

Percent  back  cavitation  allowed 

Blade  area  ratio 

Error  message  indication  as  defined 

in  Common  Block  BARCK 

P/D  of  propeller 

Internal  flag  indicating  BAR 

optimization 

f 0  if  not  optimized 

^ 1  if  optimized 


COWENTS : 


NAME: 


FUNCTION  CALC 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCK: 
SUBROUTINE  CALLED: 
CALLED  BY: 


FUNCTION  CALC  calculates  "KT"  or  "KQ"  for 
two  to  seven-bladed  propellers,  given  J, 
BAR,  P/D,  Z,  and  an  indicator  for  "KT"  or 
"KQ." 


TQ  =  CALC  (J,  BAR,  PD,  Z,  ITQ) 


J 

Advance  coefficient  of  propeller 

BAR 

Blade  area  ratio 

PD 

P/D  of  propeller 

Z 

Number  of  blades 

ITQ 

Index  defined  by 

fl  if  thrust  is  input 

ITQ  = 

1^2  if  power  is  input 

OFFDES 


NONE 


TROOST,  ROOT,  POLCOF 


COMMENTS : 


NAME: 


FUNCTION  CFITTC 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENT: 

COWON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


Function  CFITTC  determines  the  ITTC 
frictional  resistance  coefficient  Cp 

CF  =  CFITTC  (RN) 

RN  Reynolds  number  at  at  test  conditions 

NONE 

NONE 

SYNTH 

_  0.075 _ 

F  Clogi0(Rn)-2l2 


NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS : 


SUBROUTINE  CHEB 

Subroutine  CHEB  determines  the  Chebychev 
coefficients  and  wetted  surface  area  of 
the  strut  and  body.  It  also  plots  the 
strut  waterline  and  the  body  area 
distribution. 

CALL  CHEB  (TITLE) 

TITLE  Array  containing  the  alphanumeric  char¬ 
acters  of  the  title  of  the  experiment 

COEFS,  INPUT,  OMEGA,  PHYSCO 

WSURFB,  WSURFS ,  PCHEB 

SYNTH 
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I 


NAME: 

PURPOSE: 

CALLING  SEQUENCE- 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS : 


i 

: 


SUBROUTINE  ERROR 

Subroutine  ERROR  prints  the  error  message  of 
the  propulsor  design. 

CALL  ERROR  (IE) 

IE  Error  message  pointer  as  defined 
in  Common  Block  BARCK 

NONE 

NONE 

PROPP 


COMMENTS : 


Cos [(2M-1) (0)] 


NAME: 

PURPOSE: 

CALLING  SEQUENCE: 

ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS : 


SUBROUTINE  FINDRG 

Subroutine  FINDRG  determines  the  fin  drag. 

CALL  FINDRG  (RHOS,  XNUS,  VSFW,  SFINS, 
TFINS,  CFINS ,  TBFINS,  SBFINS,  DS) 

RHOS  Density  of  water 

XNUS  Kinematic  viscosity  of  water 

VSFW  Ship  speed  (fps) 

SFINS  Chord  *  span 

TFINS  Thickness  of  a  fin  (ft) 

CFINS  Chord  of  a  fin  (ft) 

TBFINS  O.CI  *  TFINS 

SBFINS  TBFINS  *  span 

DS  Total  fin  drag 

NONE 

FRICT,  FOILSC 
SYNTH 
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NAME: 

PURPOSE: 

CALLING  SEQUENCE: 

ARGUMENTS : 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 


SUBROUTINE  FOILSC 

Subroutine  FOILSC  calculates  the  drag 
components  for  a  symmetrical  foil  section 

CALL  FOILSC  (XS,  CFS,  TOC,  TOC2,  TOC4,  SS 
TS2,  SBS,  CLOTB ,  RFS,  RFVAS ,  RPS,  RINTS, 
RBS) 

XS  Ship  constant  =  h  *  p  *  SFINS 

VSFW  **  2 

CFS  Flat  plate  friction  coefficient 

TOC  Thickness/chord  ratio 

TOC2  T/C**2 

T0C4  T/C**4 

SS  Chord  *  span  -  area  of  foil 

TS2  (SS)  **  2  -  area  squared 

SBS  TBFINS  *  span 

CLOTB  Chord/0.01  *  TFINS 

RFS  Frictional  resistance  of  a  flat 

plate 

RFUAS  Resistance  due  to  velocity 

argumentation 

RPS  Resistance  due  to  pressure  and 

separation 

RINTS  Resistance  due  to  intersection 

with  hull 

RBS  Base  drag  due  to  bluntness  of 

trailing  edge 

NONE 

NONE 

FINDRG 


COMMENTS: 


NAME: 


FUNCTION  FORMDR 


PURPOSE : 


Function  FORMDR  evaluates  the  form  drag 
coefficient . 


CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


FDR  =  FORMDR  (VL) 

VL  Ratio  of  ship  speed  to  square  root  of 
length  of  strut  V//L^~ 

NONE 

FUNCTION  YINTP 
RWOUT 
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NAME: 


SUBROUTINE  FRICT 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  FRICT  calculates  the  flat  plate 
friction  coefficient  based  on  the  local 
Reynolds  number. 

CALL  FRICT  (CF,  VF,  XLEN ,  XNU) 

CF  Flat  plate  friction  coefficient 

VF  Speed  of  ship  (fps) 

XLEN  Chord  of  fin  (ft) 

XNU  Kinematic  viscosity 

NONE 

NONE 

FINDRG 


NAME: 


SUBROUTINE  PCHEB 


PURPOSE: 


CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  PCHEB  plots  by  line  printer  the 
body  sectional  area  curve  and  waterplane 
outline  curve  from  the  given  Chebychev 
coefficients. 

CALL  PCHEB  (AS,  BS,  AB,  BB,  NN,  TITLE) 

AS  Coefficients  of  Chebychev  Sine 

Series  for  strut 

BS  Coefficients  of  Chebychev  Cosine 

Series  for  strut 

AB  Coefficients  of  Chebychev  Sine 

Series  for  body 

BB  Coefficients  of  Chebychev  Cosine 

Series  for  body 

NN  Maximum  order  of  Chebychev  Series 

TITLE  Array  containing  the  alphanumeric  char¬ 

acters  of  the  title  of  the  experiment 

PLOT,  XRPLOTQ 

PLOT1 ,  PLOT2,  PLOT3,  PLOT4 
CHEB 
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NAME: 


SUBROUTINE  PLOT1 


i  l 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCK: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  PLOT1  sets  up  spacing  and 
determines  the  values  of  the  axes. 


CALL  PLOT  1  (NSCALE,  A,  B,  C,  D,  E,  F) 


NSCALE 

NSCALE  (1) 
NSCALE  (2) 
NSCALE  (3) 
NSCALE  (4) 

A 

B 

C 

D 

E 

F 

XRPLOTQ 

NONE 

PCHEB 


Integer  array  defined  as 
follows: 

I,  if  printed  values  of  the 
ordinate  are  10  **  I  times 
the  actual  value 

J,  if  printed  values  of  the 
ordinate  are  10  **  J  times 
the  actual  value 

K,  if  printed  values  of  the 
abscissa  are  10  **  K  times 
the  actual  values 

L,  if  printed  values  of  the 

abscissa  are  10  **  L  times  | 

the  actual  value 

Integer  number  of 
horizontal  grid  lines 
Integer  number  of  spaces 
beyond  each  horizontal  grid 
line  to  the  next  grid  line 
Integer  number  of  vertical 
grid  lines 

Integer  number  of  spaces 
beyond  each  vertical  grid 
line  to  the  next  grid  line 
Horizontal  grid  character 
Vertical  grid  character 


i 


NAME: 

PURPOSE: 


CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


SUBROUTINE  PLOT2 

Subroutine  PLOT2  examines  the  minimum  and 
maximum  values  of  the  abscissa  and  the 
ordinate  and  establishes  an  internal 
formula  for  computing  location  in  the  image 
region  corresponding  to  the  point  to  be 
plotted. 

CALL  PL0T2  (XMAX,  XMIN,  YMAX,  YMIN,  NSCLI) 


XMAX 

Value  of  abscissa  at  right-most 
grid  line 

XMIN 

Value  of  abscissa  at  left-most 
grid  line 

YMAX 

Value  of  ordinate  at  top  grid 
line 

YMIN 

Value  of  ordinate  at  bottom  grid 
line 

NSCLI 

Logical  flag  (should  be  false, 
if  PLOT1  has  not  been  called  and 
standard  grid  is  desired) 

XRPLOTF , 

XRPLOTQ,  XRPLOTG 

NONE 

PCHEB 
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NAME: 


SUBROUTINE  PLOT3 


PURPOSE:  Subroutine  PLOT3  assigns  an  alpha-character 

to  each  point  to  be  plotted. 

CALLING  SEQUENCE:  CALL  PL0T3  (PCHAR,  X,  Y,  SDATA,  FDATA,  DDATA) 


ARGUMENTS: 

PCHAR 

Plotting  character 

X 

Array  containing  the  X  coordinates 
to  be  plotted 

Y 

Array  containing  the  Y  coordinates 
to  be  plotted 

SDATA 

Integer  position  in  the  arrays  of 
the  first  ordered  pair  to  be  plotted 
' 1  if  each  point  from  SDATA  to  DDATA 
is  to  be  plotted 

FDATA  =< 

2  if  every  other  point  is  to  be 
plotted 

3  if  every  third  point  is  to  be 
plotted 

DDATA 

Integer  position  in  the  array  of  the 
last  ordered  pair  to  be  plotted 

COMMON  BLOCKS: 

XRPLOTF, 

XRPLOTG 

SUBROUTINE  CALLED: 

NONE 

CALLED  BY: 

COMMENTS: 

PCHEB 
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NAME: 


SUBROUTINE  PLOT4 


PURPOSE: 


CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 


Subroutine  PLOT4  prints  the  image  of  the 
completed  graph  on  the  printer,  including 
the  values  of  the  abscissa  and  the  ordinate 
at  the  grid  lines  outside  the  bottom  and 
left  edge  of  the  graph. 

CALL  PL0T4  (MCHAR,  NCHAR) 

MCHAR  Single  dimension  array  containing 

alpha  characters  to  be  plotted  at 
the  left  of  the  graph 
NCHAR  Number  of  valid  characters  in 

MCHAR 

XRPLOTF,  XRPLOTG,  XRPLOTQ 

QPLOTZ5 

PCHEB 


COMMENTS: 


NAME: 

PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


SUBROUTINE  POLCOF 

Subroutine  POLCOF  determines  the 
coefficients  of  the  polynomials  which 
approximate  the  KT  and  KQ  curves. 

CALL  POLCOF  C  PD,  BAR,  Z) 

PD  P/D  of  propeller 

BAR  Blade  area  ratio 

Z  Number  of  blades 

OFFDES,  PCOF 

Function  CALC 

PROPER 
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NAME: 


SUBROUTINE  PRODES 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS: 


Subroutine  PRODES  is  used  in  the  design 
of  a  propulsion  system  to  determine 
characteristics  and  performance  at  the 
design  condition. 

CALL  PRODES  (IS) 

IS  Type  of  propulsor  (in  this 

program  set  equal  to  1  to 
indicate  conventional  fixed 
blade  propellers) 

BARCK,  CDESC,  CPROP 

TROOST,  BURIL 

SHPCMP 


NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS: 


SUBROUTINE  PROPER 

Subroutine  PROPER  calculates  the  open  water 
curves'  table  and  off -design  performance  of 
propulsor  system  "IS." 


CALL  PROPER  (IS,  VOF,  TOF,  PWR,  TITLE) 


IS 

Type  of  propulsor  system  (in  this 
program,  set  to  1  to  indicate 
conventional  fixed  blade  propellers) 

VOF 

Off  design  speed  (ft) 

TOF 

Off  design  thrust  (lbs) 

PWR 

SHP  at  off  design  condition 

TITLE 

Array  containing  the  alphanumeric 
characters  of  the  title  of  the 

experiment 

CDESC,  CPROPR,  PCOF,  ROOTC 
POLCOF,  ROOTP 


SHPCMP 
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NAME: 

SUBROUTINE 

PROPP 

PURPOSE : 

Subroutine 

PROPP  prints  all  data  in 

common  block  "CPROP" 

CALLING  SEQUENCE: 

CALL  PROPP 

1  (IS,  TITLE) 

ARGUMENTS: 

IS 

Type  of  propul sor  system  (in  this 
program,  set  to  1  to  indicate 
conventional  fixed  blade  propellers) 

TITLE 

Array  containing  the  alphanumeric 
characters  of  the  title  of  the 
experiment 

COMMON  BLOCK: 

CPROP 

SUBROUTINE  CALLED: 

ERROR 

CALLED  BY: 

SHPCMP 

COMMENTS: 
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NAME: 


SUBROUTINE  QPLOTZ5 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
COMMENTS : 


Subroutine  QPLOTZ5  calculates  the  scaling 
information  needed  to  generate  the  format 
to  label  the  left-hand  side  of  the  program. 

CALL  QPL0TZ5  (PDQ) 

PDQ  Scaling  factor  for  ordinate  plot 

XRPLOTF ,  XRPLOTQ 

NONE 


I 


S  * 

5 

% 


NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 

COMMONG  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 


SUBROUTINE  REFLKT 

Subroutine  REFLKT  reflects  the  symmetrical 
matrices  to  T  and  W. 

CALL  REFLKT 

NONE 

AUX,  COEFS 

NONE 

RWAVE 


COMMENTS: 


NAME: 

SUBROUTINE  RIN 

PURPOSE: 

Subroutine  RIN  reads  input 
program . 

data  for  the 

CALLING  SEQUENCE: 

CALL  RIN  (TITLE) 

ARGUMENT: 

TITLE  Array  containing 

the  alphanxuneric  char 

acters  of  the  title  of  the  experiment 

COMMON  BLOCK: 

INPUT 

SUBROUTINE  CALLED: 

NONE 

CALLED  BY: 

SYNTH 

COMMENTS: 
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NAME: 

PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS: 

COfWON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


SUBROUTINE  RINIT 

Subroutine  RINIT  initializes  the  T  and  W 
arrays . 

CALL  RINIT 

NONE 

AUX,  COEFS 

NONE 

RWAVE 
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NAME: 

PURPOSE : 

CALLING  SEQUENCE: 

ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


SUBROUTINE  RINTEG 

Subroutine  RINTEG  provides  numerical 
integration  for  the  auxiliary  function 
of  T  and  W. 

CALL  RINTEG  (ALFA,  B,  D,  NLOC2 ,  WTINT, 

SEPCOS,  SQ) 

ALFA  Integrating  variable  (o) 

B  Ratio  of  distance  between 

centerlines  of  strut  and  body 
to  the  length  of  strut 
D  Ratio  of  distance  between 

centerlines  of  second  strut 
and  body  to  the  length  of 
strut 

f  0  if  single  strut 
NL0C2  =  \ 

1  if  tandem  struts 

WTINT  Weighting  constant  for  numerical 

integration 

SEPCOS  Value  of  the  Cosine  function  in 
the  integrand 

1 

SQ  Value  of  ^2_  y2)^  in  the  integrand 

AUX,  COEFS,  INPUT,  OMEGA 

BESSJ 

RWAVE 
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NAME: 


SUBROUTINE  ROOT 


PURPOSE:  Subroutine  ROOT  uses  the  Regula  Falsi 

method  to  find  the  intersection  of  the 
"K"  curve  and  the  "K/J**2"  curve. 


CALLING  SEQUENCE:  CALL  ROOT  (PD,  Z,  BAR,  COEFF,  IPJ,  ITQ, 

KT,  KQ,  EFFIN,  J2) 


ARGUMENTS: 

PD 

P/D  of  propeller 

Z 

Number  of  blades 

BAR 

Blade  area  ratio 

COEFF 

Optimization  coefficients  as 
defined  is  TROOST 

IPJ 

Optimization  index  defined  by: 
f  1 ,  no  optimization 

2,  optimize  N  with  KT/J**2 

IPJ  =  < 

3,  optimize  N  with  KQ/J**3 

4,  optimize  D  with  KT/J**4 

5,  optimize  D  with  KQ/J**5 

ITQ 

ITQ  =  < 

Index  indicating  the  following 
pi,  if  thrust  is  input 

^2,  if  power  is  input 

KT 

Thrust  coefficient 

KQ 

Torque  coefficient 

EFFIN 

Open  water  efficiency 

J2 

Advance  coefficient 

COMMON  BLOCK: 

ROOTC 

SUBROUTINE  CALLED: 

FUNCTION 

CALC 

CALLED  BY: 

TROOST 

COMMENTS : 
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NAME: 


SUBROUTINE  ROUT 


PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS : 

COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS: 


Subroutine  ROUT  outputs  the  computed  results 
for  all  components  of  drag  coefficients  and 
effective  horsepower. 

CALL  ROUT  (TITLE,  EHP) 

TITLE  Array  containing  the  alphanumeric  char¬ 
acters  of  the  title  of  the  experiment 
EHP  Effective  horsepower  of  the  ship 

INPUT,  OMEGA,  PHYSCO 

SUM,  Function  FORMDR 

SYNTH 
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NAME: 


SUBROUTINE  RWAVE 


PURPOSE:  Subroutine  RWAVE  evaluates  the  auxiliary 

functions  of  T  and  W. 

CALLING  SEQUENCE:  CALL  RWAVE  (B,  D,  NLOC2) 

ARGUMENTS:  B  Ratio  of  distance  between  center- 

lines  of  strut  and  body  to  the 
length  of  strut 

D  Ratio  of  distance  between  center- 


lines  of  second  strut  and  body  to 
the  length  of  strut 


NL0C2 

Switch  indicating  the  presence  of 
a  second  strut 

f  0  if  single  strut 

NL0C2 

*1 

1 

^ 1  if  tandem  struts 

COMMON  BLOCKS: 

AUX,  COEFS,  INPUT,  OMEGA,  PHYSCO,  PSI 

SUBROUTINES  CALLED: 

RINIT, 

SIMPSN,  RINTEG,  BESSJ,  REFLKT 

CALLED  BY: 

SWATH 

COMMENTS: 

Typical  Auxiliary  Functions 
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NAME: 


SUBROUTINE  SHPCMP 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS : 


Subroutine  SHPCMP  determines  a  propeller's 
design  and  off-design  performance  and  uses 
this  to  determine  shaft  horsepower. 

CALL  SHPCMP  (VDES,  EHPDES,  PDIA,  BDIA,  HB, 

VOFF,  EHP,  SHP,  NV,  TITLE) 

VDES  Ship  speed  at  design  point  (fps) 

EHPDES  Effective  horsepower  at  design 

point 

PDIA  Diameter  of  propeller  (ft) 

BDIA  Diameter  of  hull  body  (ft) 

HB  Draft  of  body  (ft) 

VOFF  Off-design  speed  of  ship  (fps) 

EHP  Effect ‘ ve  horsepower  at  off 

«*osign  speeds 

SHP  Shaft  horsepower  at  off  design 

speeds 

NV  Numbvr  of  test  speeds 

TITLE  Array  containing  the  alphanumeric  char- 

.?.*  ters  of  the  title  of  the  experiment 

CDESC,  CPROP 

PRODES ,  PROPP ,  PROPER,  PROPERO 
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NAME: 


SUBROUTINE  SIMPSN 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


Subroutine  SIMPSN  sets  up  Simpson's 
multiplier  table  for  numerical  integration 
using  Simpson's  rule. 

CALL  SIMPSN  (NPTS,  SIMP) 


NPTS 

Number  of  integration  steps  (odd 
and  greater  than  one) 

SIMP 

Array  containing  Simpson's 
coefficients 

NONE 

NONE 

RWAVE 
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NAME: 


SUBROUTINE  SUM 


PURPOSE: 


CALLING  SEQUENCE: 
ARGUMENTS : 


COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


Subroutine  SUM  computes  the  matrix  products 

as  follows: 

n  m 

SUM  I  I  {Asin  Asn  Tgmn  +  Bsm  Bsn  Wsmn) 

m=l  n=l 

CALL  SUM  (SUM1S,  SUM1B,  SUM1SB,  SUM12S , 

SUM12B,  SUM125B) 

SUM1S  Partial  sum  for  strut  1 

SUM1B  Partial  sum  for  body  1 

SUM1SB  Partial  sum  for  interaction  between 

strut  1  and  body  1 

SUBM12S  Partial  sum  for  interaction  between 

strut  1  and  strut  2 

SUM12B  Partial  sum  for  interaction  between 

body  1  and  body  2 

SUM12SB  Partial  sum  for  interactions  between 
body  1  and  strut  2  or  body  2  and 
strut  1 


AUX,  COEFS 

NONE 

ROUT 
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NAME: 

PURPOSE : 

CALLING  SEQUENCE: 

ARGUMENTS: 


COMMON  BLOCKS: 
SUBROUTINES  CALLED: 
CALLED  BY: 

COMMENTS : 


SUBROUTINE  TROOST 

Subroutine  TROOST  is  used  with  Subroutine 
ROOT  and  Function  CALC.  Either  rpm  or 
diameter  can  be  optimized  given  the  other. 

CALL  TROOST  (VAI,  P,  T,  Nl,  D,  Z,  PDNEW , 
BAR,  KT,  KQ,  EFFNEW,  J,  ITQ,  IPJ,  RHO, 
IERROR) 


VAI 

Waked  ship  velocity  (kts) 

P 

Torque  per  shaft  (ft- lb) 

T 

Propeller  thrust  (lbs  per  shaft) 

Nl 

Propeller  rpm 

D 

Propeller  diameter  (ft) 

Z 

Number  of  blades  of  propeller 

PDNEW 

P/D  of  propeller 

BAR 

Blade  area  ratio 

KT 

Thrust  coefficient 

KQ 

Torque  coefficient 

EFFNEW 

Open  water  efficiency 

J 

Advance  coefficient 

ITQ 

Index  defined  by: 

ITQ  =  s 

IPJ 

1  if  thrust  is  input 

2  if  power  is  input 

Optimization  index  defined  by: 

j 

1  -  no  optimization,  J  known. 

i 

IPJ  =  \ 

calculates  KT,  KQ,  EFFNEW  on 
2  -  optimize  N  with  KT/J**2 
,  3  -  optimize  N  with  KQ/J**3 

RHO 

4  -  optimize  D  with  KT/J**4 
v  5  -  optimize  D  with  KQ/J**5 

Water  density 

IERROR 

Error  message  indicator 

NONE 


ROOT,  TR0ST2 ,  FUNCTION  CALC 
PRODES 
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NAME: 


SUBROUTINE  TROST2 


PURPOSE : 

CALLING  SEQUENCE: 

ARGUMENTS : 


COKMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 
COMMENTS: 


Subroutine  TROST2  is  used  to  calculate 
a  single  KT,  KQ  and  efficiency  value  for 
a  given  propeller. 

CALL  TR0ST2  (KT,  KQ,  J,  PDTR2 ,  BAR,  Z, 

EFFTR2) 

KT  Thrust  coefficient 

KQ  Torque  coefficient 

J  Advance  coefficient 

PDTR2  Pitch  to  diameter  ratio 

BAR  Blade  area  ratio 

Z  Number  of  blades 

EFFTR2  Open  water  efficiency 

NONE  1 

Function  CALC 

TROOST  I 
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NAME: 


SUBROUTINE  WSURFB 


PURPOSE:  Subroutine  WSURFB  determines  the  surface 

area  of  a  body  of  revolution  given  by 
Chebychev  coefficients. 

CALLING  SEQUENCE:  CALL  WFSURFB  (AREA) 

ARGUMENT:  AREA  Surface  area  of  a  body  (ft^) 

COMMON  BLOCKS:  COEFS,  INPUT,  PHYSCO 

SUBROUTINE  CALLED:  FUNCTION  EVAL 

CALLED  BY:  CHEB 

COMMENTS : 
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NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENT: 

COMMON  BLOCKS: 
SUBROUTINE  CALLED: 
CALLED  BY: 

COMMENTS : 


V 

I 


SUBROUTINE  WSURFS 

Subroutine  WSURFS  determines  the  wetted 
surface  area  of  a  strut. 

CALL  WSURFS  (AREA) 

AREA  Wetted  surface  area  of  a  strut 
COEFS,  INPUT 
FUNCTION  EVAL 
CHEB 


*  + 


NAME: 


FUNCTION  YINTP 


PURPOSE: 

Function 

YINTP  interpolates  through 

of  discrete  data. 

CALLING  SEQUENCE: 

YTP  = YINTP  (XA,  X,  Y,  N) 

ARGUMENTS: 

XA 

Point  to  be  interpolated 

X 

Array  of  ordinate  data 

Y 

Array  of  abscissa  data 

N 

Number  of  data  points 

COMMON  BLOCKS: 

NONE 

SUBROUTINE  CALLED: 

NONE 

CALLED  BY: 

FORMDR 

COMMENTS : 

REFERENCES 


(1)  Lin,  L.W.  and  W.G  Day,  "The  Still-Water  Resistance  and  Propulsion 
Characteristics  of  Small-Waterplane-Area  Twin-Hull  (SWATH)  Ships." 
AIAA/SNAME  Advanced  Marine  Vehicles  Conference.  San  Diego, 
California,  1974. 

(2)  Proceedings  of  the  10th  International  Towing  Tank  Conference, 
London,  England,  1963.  Published  by  the  National  Physical 
Laboratory,  England. 

(3)  Lasky,  M.P.,  "An  Investigation  of  Appendage  Drag.”  NSRDC/SPD 
Report  458-H-01 . 

(4)  Roddy,  R.F.  and  J.  Strom-Tej sen ,  "Study  of  Air  Cushion  Vehicle 
Propulsion  System  -  Propulsor  Design  and  Performance."  NSRDC/SPD 
Report  378-09. 
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PROGRAM  LISTING  OF  CLOSEFIT 
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BLOCK  DATA  BLDATA 

BLDATA 

C0MM0N/0UT/HS,HB,XLS.XL3,TSMAX,AX,PI , G, RHO. GNU , WETS, WETB.WTSURF,  BLDATA 
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DTNSROC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSROC  REPORTS.  A  FORMAL  SERIES.  CONTAIN  INFORMATION  OF  PERMANENT  TECH 

NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF  \ 

THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM 
INARY,  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 

THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INIEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSROC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


